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PHASE TRANSITIONS IN DELAUNAY POTTS MODELS 


STEFAN ADAMS AND MICHAEL EYERS 


Abstract. We establish phase transitions for certain classes of continuum Delaunay multi¬ 
type particle systems (continuum Potts models) with infinite range repulsive interaction be¬ 
tween particles of different type. In one class of the Delaunay Potts models studied the repulsive 
interaction is a triangle (multi-body) interaction whereas in the second class the interaction is 
between pairs (edges) of the Delaunay graph. The result for the edge model is an extension 
of finite range results in 1BBD04] for the Delaunay graph and in GH96] for continuum Potts 
models to an infinite range repulsion decaying with the edge length. This is a proof of an 
old conjecture of Lebowitz and Lieb. The repulsive triangle interactions have infinite range 
as well and depend on the underlying geometry and thus are a first step towards studying 
phase transitions for geometry-dependent multi-body systems. Our approach involves a De¬ 
launay random-cluster representation analogous to the Fortuin-Kasteleyn representation of the 
Potts model. The phase transitions manifest themselves in the percolation of the correspond¬ 
ing random-cluster model. Our proofs rely on recent studies DDG12] of Gibbs measures for 
geometry-dependent interactions. 


1. Introduction 

Although the study of phase transitions is one of the main subjects of mathematical statis¬ 
tical mechanics, examples of models exhibiting phase transition are mainly restricted to lattice 
systems. In the continuous setting results are much harder to obtain, e.g., the proof of a liquid- 
vapor phase transition in LMP99] , or the spontaneous breaking of rotational symmetry in two 
dimensions for a Delaunay hard-equilaterality like interaction [ MR09] . These phase transitions 
manifest themselves in breaking of a continuous symmetry. There is another specific model 
for which a phase transition is known to occur: the model of Widom and Rowlinson [ WR70 ]. 
This is a multi-type particle system in M d , d > 2, with hard-core exclusion between particles 
of different type, and no interaction between particles of the same type. The phase transition 
in this model was stablished by Ruellc [Rue m Lebowitz and Lieb [ LL72 J extended his result 
by replacing the hard-core exclusion by a soft-core repulsion between unlike particles. Finally, 
phase transition results for a general class of continuum Potts models in R d , d > 2, have been 
derived in |GH96] . The repulsive interaction between particles of different type in [ GH96 ] is of 
finite range and a type-independent background potential has been added. The phase transi¬ 
tions for large activities in all these systems reveal themselves in breaking of the symmetry in 
the type-distribution. The results in [ GH96] are based an a random-cluster representation anal¬ 
ogous to the Fortuin-Kasteleyn representation of lattice Potts models, see [GHM] , In |BBD04| 
the soft repulsion in |GH96 j between unlike particles has been replaced by another kind of 
soft repulsion based on the structure of some graph. More precisely, the finite range repulsion 
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between particles of different type acts only on a nearest neighbour subgraph of the Delaunay 
graph in M 2 . 


In this paper we establish the existence of phase transitions for two different classes of con¬ 
tinuum Delaunay Potts models in M 2 including the Widom-Rowlinson model. The repulsive 
interaction between unlike particles in all models is of infinite range, and it depends in one class 
on the geometry of the Delaunay triangle (minimal angle) and in another class on the Delau¬ 
nay edge when the interaction decays with the length of the edge. All our models come with a 
hard-core background interaction which is independent of the type of the particles. We believe 
that one can extend our results by removing this hard-core constraint, see Remark 2/T below. 
Our results extend [ BBD04] and [LGH96] in two ways by having an infinite range repulsion and a 
geometry-dependency of the type interaction. Gibbs models on Delaunay structures have been 


studied in [BBD991 IBBD021 IBBD041 IDerOSl IDDG121 IDG091 IDLllj . and our results rely on 
the existence of Gibbs measures for geometry-dependent interactions established in [ DDG12] . 
Our approach is based on a Delaunay random-cluster representation similar to |BBD04] and 
[GH96j . the difference being that we replace edge percolation by an adaptation of lattice hyper¬ 
edge percolation )Gri94j to our continuum setting. A phase transition for our Delaunay Potts 
models follows if we can show that the corresponding percolation process contains an infinite 
cluster. A similar program was carried out by Chayes et al. in [CCK95] for the hard-core 
Widom-Rowlinson model. In that case, the existence of infinite clusters follows from a stochas¬ 
tic comparison with the Poisson Boolean model of continuum percolation, while our framework 
uses a coarse graining method to derive a stochastic comparison with site percolation on 1?. 
This idea actually goes back to |HagOO| in which percolation for the Poisson Voronoi model has 
been established. We note that our random-cluster representation requires the symmetry of 
the type interaction. In the non-symmetric Widom-Rowlinson models, the existence of a phase 
transition has been established by Bricmont et al. (BKL] . and recently by Suhov et al. [ MSS ]. 


We conclude with some remarks on the particular features models defined on the Delaunay 
hypergraph structure show. The most simple Delaunay Potts model would be one with no 
more interaction than constant Delaunay edge interaction as then the percolation would be 
independent of the activity parameter. Our triangle model is a step towards understanding 
this system. We draw attention to some differences between geometric models on the Delau¬ 
nay hypergraph structure and that of classical models such as the Widom-Rowlinson model 
and its soft-core variant of Lebowitz and Lieb jLL72] . The first is that edges and triangles 
in the Delaunay hypergraph are each proportional in number to the number of particles in 
the configuration. However, in the case of the complete hypergraph the number of edges is 
proportional to the number of particles squared and the number of triangles is proportional 
to the number of particles cubed. Secondly, in the complete graph of the classical models, 
the neighbourhood of a given point depends only on the distance between points and so the 
number of neighbours increases with the activity parameter z of the underlying point process. 
This means that the system will become strongly connected for high values of z. This is not 
the case for the Delaunay hypergraphs which exhibit a self-similar property. Essentially, as 
the activity parameter z increases, the expected number of neighbours to a given point in the 
Delaunay hypergraph remains the same, see |Mp94| . Therefore, in order to keep a strong con¬ 
nectivity in our geometric models on Delaunay hypergraphs, we use a type interaction between 
particles of a hyperedge with a non-constant mark. Finally, and perhaps most importantly, is 
the question of additivity. Namely, suppose we have an existing particle configuration u and 
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we want to add a new particle x to it. In the case of classical many-body interactions, this 
addition will introduce new interactions that occur between x and the existing configuration oj. 
However, the interactions between particles of oj remain unaffected, and so classical many-body 
interactions are additive. On the other hand, in the Delaunay framework, the introduction of a 
new particle to an existing configuration not only creates new edges and triangles, but destroys 
some too. The Delaunay interactions are therefore not additive, and for this reason, attractive 
and repulsive interactions are indistinct. In the case of a hard exclusion interaction, we arrive 
at the possibility that a configuration oj is excluded, but for some x, oj U x is not. This is 
called the non-hereditary property |DG09j . which seems to rule out using techniques such as 
stochastic comparisons of point processes GK97j . 

2. Results 

2.1. Setup. We consider configurations of points in R 2 with internal degrees of freedom, or 
marks. Let E = {1,..., q},q G N, q > 2, be the finite set of different marks. That is, each 
marked point is represented by a position x G R 2 and a mark cr(x) G E , and each marked 
configuration is a countable subset of R 2 x E having a locally finite projection onto R 2 . We 
denote by 12 the set of all marked configurations with locally finite projection onto R 2 . We will 
sometimes identify oj with a vector oj = (o; (1) , ... ,oj m ) of pairwise disjoint locally finite sets 
cu (1) ,..., a; t9) in R 2 (we write 12 for the set of all locally finite configurations in R 2 ). Any oj is 
uniquely determined by the pair (cu,c), where oj = U^ =1 a; (i, is the set of all occupied positions, 
and where the mark function a: oj —» E is defined by cr(x) = i if t G cu (i) ,i G E. For each 
measurable set B in R 2 x E the counting variable N(B ): oj —>■ w{B) on 12 gives the number of 
marked particles such that the pair (position, mark) belongs to B. We equip the space 12 with 
the cr-algebra J- generated by the counting variables N(B) and the space 12 of locally finite 
configurations with the cr-algebra T generated by the counting variables Na = fl A} for 
A d R 2 where we write A d R 2 for any bounded A C R 2 . As usual, we take as reference 
measure on (12, J-) the marked Poisson point process II 1 with intensity measure zLeb (8) p u 
where z > 0 is an arbitrary activity, Leb is the Lebesgue measure in R 2 , and /i u is the uniform 
probability measure on E. 

We let 12/ C 12 (resp. 12/ C 12) denote the set of all finite configurations. For each A C R 2 we 
write 12 a = {u;Gl2:u;cAx E} for the set of configurations in A, pr A : oj —> oja := oj fl A x E 
for the projection from 12 to 12 a (similarly for unmarked configurations), J-' K = J-\ n A for the 
trace cr-algebra of T on 12a, and Ta = C T for the a -algebra of all events that happen 

in A only. The reference measure on (12a, ) is U Z A := ITopr/. In a similar way we define the 

corresponding objects for unmarked configurations, n 2 , n^, 12a, prA; *^a, and Ea. Finally, let 
O = ('dx)xeR 2 be fb e shift group, where d x \ 12 —^ 12 is the translation of the spatial component 
by the vector — x G R 2 . Note that by definition, AA^ztj) = Na+ x (u) for all A C R 2 . 

The interaction between the points in all the models to be studied depend on the geometry 
of their location. We describe this in terms of hypergraph structures. A hypergraph structure 
is a measurable set of 12/ x 12 (resp. 12/ x 12). We outline the definitions for the unmarked 
configurations first with obvious adaptations to the case of marked point configurations. The 
set Del of Delaunay hyperedges consist of all pairs with rj C oj for which there exits 

an open ball B(rj,oj) with dB(rj,oj) oj = r) that contains no points of oj. For m = 1,2,3, 
we write Del m = G Del: = m} for the set of Delaunay simplices with m vertices. 
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Given a configuration oj the set of all Delaunay hyperedges y C oj with #r] = m is denoted 
by Del m (u;). It is possible that r/ G Del(cu) consists of four or more points on a sphere with 
no points inside. In fact, for this not to happen, we must consider configurations in general 
position as in |M094| . More precisely, this means that no four points lie on the boundary of a 
circle and every half-plane contains at least one point. Fortunately, this occurs with probability 
one for our Poisson reference measures, and in fact, for any stationary point process. Note that 
the open ball B(t),oj ) is only uniquely determined when jpy = 3 and rj is affinely independent. 
Henceforth, for each configuration oj we have an associated Delaunay triangulation 

{r C oj: = 3, B(t, oj) fl oj = 0} (2.1) 

of the plane, where B(t,oj) is the unique open ball with r C dB(r,ui). The set in ( |2.1 ) 
is uniquely determined and defines a triangulation of the convex hull of oj whenever oj is in 
general position ( |M094| ). In a similar way one can define the marked Delaunay hyperedges, 
i.e., Del and Del m (cu) as measurable sets in ft f x it where the Delaunay property refers to the 
spatial component only. 


Given a configuration u 6 0 (or oj) we write = {<CgD:<C\A = ci;} (resp. fi\A for 
the set of configurations which equal oj off A. For any triangle r G Del 3 we denote the minimal 
angle by a(r) G (0, 7 t/ 3 ] and for any edge y G Del 2 we denote its length by £(y). We write 
Adi 2 in the following when A is non-empty and finite, i.e., 0 < |A| < 00 . The interaction is 
given by the following Hamiltonian in A with boundary condition oj G ft, m = 2, 3, 

HaAC) : = x-y)+ ^ (p{a(r))(l-5 a {r)), C e ^a, u , (2.2) 

7?EDel2 a(C) ' r ^ Del m(C)! 

rGDel mA (C) 


where Del 2 ,A(C) : = \ r l G Del 2 ((): 3r G Del 3 (C), V C r, 8B(t, () fl A ^ 0 } and Del 3jA (C) : = ( r e 
Del 3 (<C) : 0B(t, () fl A 7 ^ 0 }. Here ip: M 2 —>• M U {c>o} is an even measurable function and 0 is 
a measurable function of the minimal angle a(r) when m — 3 or a measurable function of the 
length £(rj) of an edge when m = 2 (see below) and 


SJt) = 


1 , if a T (x) = cr T (y) for any pair {x,y} C r, 

0 , otherwise. 


The second term on the right hand side of (2.2) is the type interaction and thus the most 
important contribution. It describes a repulsion between particles having a different mark. 
The first term corresponds to a ’colour-blind’ background interaction. Specifically, we assume 
the following for our models. 


Assumptions: 

(Al) (i) Triangle-model I: m — 3 and 0: (0, tt/3] —> M + and there are «o £ (0, 7 r/ 3 ) and 
/3 > 0 such that 

0 , if 9 < a 0 , 

/3 , if 9 > a 0 - 

(ii) Triangle-model II: m — 3 and 0: (0, 7 t/ 3 ] —> M + with 

0(0)=log(l + 00 3 ), 0> 0. 


m = 
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(iii) Edge-model: Let <5o > 0 be some parameter to be specified in (A2), m = 2 and 
0: [0, oo] —> M+ U {cxo} with 

0(^) = log(l + f3(5 0 /£) 3 ), f3 > 0. 

(A2) -0 is a hard-core potential, that is, there is a range 5 0 > 0 such that if(x — y) = 0 
whenever \x — y\ > <5 0 and ^(x — y) = +oo otherwise. 


Following [ DDtjT2j we say a configuration lo e f2 (or tu G hi) is admissible for A d M 2 and 
activity z if H A ^ is II''-almost surely well-defined and 0 < Z A (co) < oo, where the partition 
function is defined as 


ZaH 



e -^A,Ao n 2 (dc (1 ) )---n 2 (dc w ). 


We denote the set of admissible configurations by fl* A . The Gibbs distribution for y ;, 0 , and 
z > 0 in A with admissible boundary condition uj is defined as 


7a, u, (A) = 


Za(u\ 


[ ll A (CUw)e~ HA -“ (<:) n 2 (dC), AeJF. 

j ^ A UJ 


(2.3) 


It is evident from (2.3) that, for fixed ( 
C = (C (1) , • • ■, C (q) ) relative to 7 a,. 


€ the conditional distribution of the marks of 
is that of a discrete Potts model on ( embedded in the 


Delaunay triangulation with posit ion-dependent interaction between the marks. This justifies 
calling our models Delaunay Potts models. 


Definition 2.1. A probability measure p on fl is called a Gibbs measure for the Delaunay Potts 
model with activity z > 0 and interaction potentials 0 and 0 if p(fl* A ) = 1 and 

E„[/]= f jf- [ /(CUuOe-^-lOnKdOMdu,) (2.4) 

for every A d M 2 and every measurable function f. 


The equations in (2.4) are the DLR equations (after Dobrushin, Lanford, and Ruclle). 
They express that the Gibbs distribution in (2.3) is a version of the conditional probability 
p(A\3FA c ){y)). The measurability of all objects is established in |Eye 14. !DDG12 |. 


2.2. Results and remarks. We outline the results for our triangle and edge models all of 
which are generic examples of a more general class of coarse grain ready potentials, see Section[4] 
for further details. 

Proposition 2.2 (Existence of Gibbs measures). There exist at least one Gibbs measure 
for the following Delaunay Potts models with hard-core background potential of range ho > 0. 

(a) (i) Triangle model I: For any z > 0 and for any «o 7 (0,7r/3) and (3 > 0. 

(ii) Triangle model II: For any z > 0 and (3 > 0. 

(b) Edge model: For any z > 0 and (3 > 0. 

Remark 2.3 (Gibbs measures). We obtain the existence of Gibbs measures for large enough 
activity only. This is mainly due to the techniques used in jDDG12j for the existence of Gibbs 
measures h[ Geo88j ). The restriction originates from the need to obtain an upper bound for the 
Hamiltonian for some chosen configurations, namely, the so-called pseudo-periodic configura¬ 
tions (see Appendix [A| or [DDG12] ). Large activities then ensure that these pseudo-periodic 
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configurations have sufficient mass under the reference process. Existence of Gibbs measures 
for related Delaunay models have been obtained in [ BBD99. Der08l IDG09] . Note that for q = 1 
our models have no marks and Gibbs measures do exist as well (| DDG12 ]). 


Gibbs measures for the Delaunay Potts model do exist for large enough activity z > 0. A 
phase transition is said to occur if there exists more than one Gibbs measure for the Delaunay 
Potts model. The following two theorems show that this happens when the activity z and the 
parameter f3 are sufficiently large. Note that ft is a parameter for the type interaction and not 
the usual inverse temperature. The first theorem extends the results in |BBD04] for finite range 
type interaction to the case of infinite range type interaction. 

Theorem 2.4 (Phase transition - Edge model). For all 5 0 > 0 there is (3 0 = fio(^o) sueh 
that for all f3 > there is zo = zq{(3, So) such that for all z > zq, there exit at least q different 
Gibbs measures for the Delaunay edge model. 


Remark 2.5 (Range of type-interaction). The edge model is an infinite range type- 
interaction while keeping a hard-core background potential as in [BBD04] . where it is used 
on one hand to obtain a bound on the number of connected components in the hypergraph Del 2 
once a point is inserted into a configuration, and on the other hand, to bound the number of 
points in a coarse graining scheme. We have replaced the bound on the number of connected 
components by a bound of the expectation under the edge drawing measure using the hard-core 
background potential and an sufficient decay of the type interaction. We believe that one can 
relax the hard-core assumption, see Remark\2ffi. 


The following theorem shows that phase transitions occur in our Delaunay triangle models. 
As outlined in [ DG091 Remark 3.11] a triangle model similar to the edge model in [ BBD04] 
shows phase transition as well. The triangle interaction in (DG09 j is of finite range as well. 
Our triangle models have an infinite range type interaction. We obtain phase transitions for 
sufficiently large activity and parameter f3 for two models with infinite range type interactions. 
The idea for the triangle model I is that tiles (triangles) whose minimal angles exceed a given 
threshold «o gain in probability when their vertices are of the same type. The energy of a 
triangle with minimal angle is given by the constant value of the type potential </>. In our 
second triangle model II, however, this incentive is not constant as it grows with the minimal 
angle and it decays with shrinking minimal angle. 

Theorem 2.6 (Phase transition - Triangle model). The following holds for the triangle 
models. 


(a) Triangle Model I: For all 4o > 0 and «o sufficiently small there are zo = zq(S, a o) and 
(3 0 = fto(8o, a o) > 0 such that for all z > z 0 and (3 > /3 0 there exists at least q different 
Gibbs measures for the Delaunay triangle model I. 

(b) Triangle Model II: For all 4o > 0 there is = fio(ho) such that for all (3 > fio there 
is z 0 = z 0 (/3,S o) such that for all z > zq, there exit at least q different Gibbs measures 
for the Delaunay triangle model II. 


For the triangle model I note that our proof techniques require the minimal angle threshold 

In contrast to this, the triangle model II allows 


«o to be sufficiently small, see Lemma 5.1 


the smallest angles to take any value as smaller values experience vanishing preference for type 
alignment in a triangle. 
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Remark 2.7 (Hardcore versus softcore background potential). In all our models we add 
a hard-core background potential between pairs of particles. This is for technical reasons only, 
and we believe that it is possible to replace it by some triangle potential excluding minimal angles 
below a given threshold (see mm), or by the hard-equilaterality model / l)l)G 12' . (DLll j ). 
These triangle background potentials are of infinite range in terms of Delaunay edges (pairs of 
particles). They allow to bound easily the number of connected components in the Delaunay 
random-cluster model once a new point is inserted. On the other hand inserting a point in an 
empty cell now depends on all the points in neighbouring cells. An adaptation of our coarse 
graining method in Section and restricting to certain admissible boundary conditions for the 
cells allows to estimate the probability of injecting a single particle. Another possibility is to 
replace the hard-core background pair potential by a superstable pair potential (soft-core back¬ 
ground) as in |Rue70| . |GH96j . and [ BBD02 J. This allows to control the number of particles in 
cells and to estimate the probability of insertion of single points, see e.g. [CTFTMj . However, 
for the final percolation argument it seems an additional multi scale approach is needed, which 
is based in our coarse graining technique in Section |^} Suppose that in each of the sub cells 
(paralleloptopes) there is a region with a finite number of points. Then, due to the chosen 
partition ofM. 2 , there will be Delaunay triangles connecting these regions in the different cells. 
These triangles have then a minimal angle and will thus be open with high probability. Thus all 
the regions in the different cells are connected with open triangles leaving the questions whether 
within these regions all triangle are open as well. We think that this can be achieved by some 
kind a renormalisation approach with a combined bond-site percolation as in (GH96j . We hope 
to approach this idea using Delaunay potentials studied in |BBG12] replacing the standard pair 
potential in the future ( ]A15j ). 

Remark 2.8 (Free energy). One may wonder if the phase transitions manifest itself thermo¬ 
dynamically by a non-differentiability (’’discontinuity”) of the free energy (pressure). Using the 
techniques from [ Geo94j and [ DG09| . it should be possible to obtain a variational representation 
of the free energy, see also [ ACKllj for free energy representations for marked configurations. 
Then a discontinuity of the free energy can be established using our results above. For continuum 
Potts models this has been established in [ GH96 , Remark 4.3]. For a class of bounded trian¬ 
gle potential [DG09] shows that the Gibbs measures are minimisers of the free energy. The free 
energy is the level-3 large deviation rate functional for empirical processes, see |Geo94l IXHKTT1 . 


Remark 2.9 (Uniqueness of Gibbs measures). To establish uniqueness of the Gibbs mea¬ 
sure in our Delaunay Potts models one can use the Delaunay random-cluster measure C\ niU} , to 
be defined in (3.3) below. In [ GEM , Theorem 6.10] uniqueness is established once the probability 
of an open connection of the origin to infinity is vanishing for the limiting lattice version of the 
random-cluster measure, that is, for some set A d M 2 containing the origin, 


Bui C' An ,o;(A <—» A c n ) = 0, 

n—>oo 

for a sequence of boxes A n d M 2 with A n f M 2 as n —)■ oo. One way to achieve this, is 
to obtain an stochastic domination of the Delaunay random-cluster measure by the so-called 
random Delaunay edge model of hard-core particles. Using jBBB02] we know that the critical 
probabilities for both, the site and bond percolation on the Delaunay graph, are bounded from 
below. Extension to our tile (hyperedge) percolation using [Gri94j can provide a corresponding 
lower bound as well. Thus, if our parameter j3 is chosen sufficiently small, then there is no 
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percolation in our Delaunay random-cluster measures and therefore uniqueness of the Gibbs 
measure. 

It goes without saying that the results addressed here are merely a first step towards a 
closer study of phase transitions with geometry-dependent interactions. The study for Widom- 
Rowlinson or Potts models with geometry-dependent interaction is by far not complete, one 
may wish to extend the single tile (edge or triangle) interaction to mutual adjacent Voronoi cell 
interaction. The common feature of all these ’ferromagnetic’ systems is that phase transitions 
are due to breaking the symmetry of the type distribution. Breaking of continuous symmetries 
is a much harder business, see e.g. |MR09] . which shows a breaking of rotational symmetry in 
two dimensions, and which can be seen as a model of oriented particles in M 2 with a Delaunay 
hard-equilaterality interaction. Models studied in [ IDDG12 ] may be natural candidates for the 
existence of a crystallisation transition. A recent ground state study (zero temperature) [BPT| 
for multi-body interactions of Wasserstein types shows optimality of the triangular lattice. It 
seems to be promising to analyse this model for non-zero temperature in terms of crystallisation 
transition. We will investigate this further. 

Another direction for the class of ’ferromagnetic’ models with geometry-dependent interac¬ 
tion is to analyse closer the percolation phenomena, and in particular to study the conformal 
invariance of crossing probabilities as done for the Voronoi percolation (| BS98 . BR06| . and 

nm 

Last but not least it is interesting to study the higher dimensional cases in M d , d > 3, as well 
(see [ GH96 ]). We believe that our results still hold for these cases but there are some technical 
issues related to the general quadratic position such that we defer that analysis for future study. 

The rest of the paper is organised a follows. In Section [3] we define the Delaunay random- 
cluster measure for tiles (edges/triangles), and in Section [4] we prove tile percolation for our 
different Delaunay random-cluster measures for sufficiently large activities and sufficiently large 
/3. In Section [5] we gives details of our proofs. 

3. Delaunay Random Cluster measure 

For A <s M 2 and parameters z,if,(t>, and f3 we define a joint distribution of the Delaunay 
Potts model and a tile (edge/triangle) process which we call Delaunay random-cluster model. 
Here tile is used as a general name for any hyperedge of a Delaunay graph, and it will be clear 
from the context if we refer to Delaunay edges or triangles. We follow ideas in [ GH96] and 
will replace edge percolation by tile (i.e., triangle) percolation. The basic idea is to introduce 
random tiles between points in the plane. Let 

Tr 2 = {t = {i,!/,4cI 2 :i/!/,^2,z/!/} 

be the set of all possible tiles of triples of points in R 2 , likewise, let T\ be the set of all tiles 
in A and 7} for the set of tiles in ( G £Ia,u- We identify u> with (n (1) and co = (w (1) ,0,...,0). 
This allows only monochromatic boundary conditions whereas the general version involves the 
so-called Edwards-Sokal coupling (see [ GHM J for lattice Potts models). We restrict ourself to 
the former case for ease of notation. We write 

T={Tc T r 2 : T locally finite} 
for the set of all locally finite tile configurations. 
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The joint distribution is built from the following three components. 

The point distribution is given by the Gibbs distribution P z Auj in A with admissible boundary 
condition uj G kl A , interaction - 0 , and activity zq, i.e., 


o = 


Z A (u) 


exp 


E V>(r(*j)))ny(dO, 

77GDel2(C) 


(3.1) 


with ( G n A>U] and where £(rj) = \x — y\ is the length of the Delaunay edge rj = {x,y}. Note 
that for q = 1 this measure coincides with (2.3). 


The type picking mechanism for a given configuration £ G is the distribution A^a of 

the mark vector a G EA Here (cr(x)) xe £ are independent and uniformly distributed random 
variables on E with a(x) = 1 for all x G Ca c = oj. The latter condition ensures that all points 
outside of A carry the given fixed mark. 


The tile drawing mechanism. Given a point configuration Q G Da,uj we let y^ A be the distribution 
of the random tile configuration {r G T^: v(r) = 1} G 7 with the hyperedge configuration 
v G {0, l} Tc having probability 

n p(r)" M ( 1 - pItT-M 

r S T c 


with 


p(r) 


1) 


(1 - e ^ (T) )]l De i 3 ( C )(T) if r G T r2 \ T A ?, 
ilDei 3 (C)(' r ) if t G T A c. 


(3.2) 


The measure y^ A is a point process on T r2 . Note that ( —> \^ A and ( —>■ y^ A are probability 
kernels (see |Eyel4| ). Let the measure 

f>*«„(dC,dT) = -l_/-^(dOA CA (dOp c ,A(dT) 

Zj A [u) 

be supported on the set of all (£, T ) with 0 G A2a,u> and T C T0 We shall condition on the 
event that the marks of the points are constant on each connected component in the hypergraph 
(C,T fl T(0. Two distinct vertices x and y in the hypergraph are adjacent to one another if 
there exists rGf( such that {x, y} C r. A path in the hypergraph (0, T fl 70) is an alternating 
sequence vi,ti,v 2 ,t 2 , ■ ■ ■ of distinct vertices Vi and hyperedges (tiles) t 3 such that {vi, v i+ i } C t t 
for all i > 1. We write 

V{(CT)eflxI: £(1 - $„(t)) = 0} 

tGT 

for the set of marked point configurations such that all vertices of the tiles carry the same mark. 
The set A is measurable which one can see from writing the condition in the following way 
q 

5Z (V^X 1 - %>(£)) + \d) (s/)(i - %o(X)) + 31 C (o (a^)(l - il c « (j/))) = 0 

r ={ x , y , z}ET i =1 

and using the fact that (0,x) t-A 10(0 (x),i = 1,... ,q, are measurable (see [ GH961 Chapter 2]). 
Furthermore, P AjU (A) > 0, which follows easily observing P a ^(A) > P Aoj {{cu}) = e~ zq ^/Z a (oj) 
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where uj is the configuration which equals u outside of A and which is empty inside A. Hence¬ 
forth, the random-cluster measure 


p = nu 


D) 


is well-defined. As in [ GH96] we obtain the following two measures from the random-cluster 
measure P, namely if we disregard the tiles we obtain the Delaunay Gibbs distribution 7 a, a, 
in (2.3) (see |Eyel4| ). For the second measure consider the mapping sp: (£, T) —> (£, T) from 
12 x 7 onto SlxT where £ f-A ( = U^ =1 ( (i> . For each (£, T ) with T C we let K((, T ) denote 
the number of connected components in the hypergraph (£, T). The Delaunay random-cluster 
distribution on 0 x 7 is defined by 


C A ,„(dC,dT) = —^ci^P^dO^AdT), 


where P Au] is (3.1) with activity z replacing zq and where 


(3.3) 


ZaM = / [ q KIX - T, PU dC)^c.A(dT) 

is the normalisation. It is straightforward to show that P o sp -1 = (see |Eye l4|). To 
relate the influence of the boundary condition on the mark of a single point to the connectivity 
probabilities in the random-cluster model we follow [ GH96j . For any A C A, s G E, (, E 
and (C, T), with T C T c we define 


JVa,,(0 = #{C W n A}. 


Then 

AA«a c (C? T) = #{x G C n A: x belongs to a cluster connected to A c 
in T fl Del 3 (C)} 

is the number of points in £ fl A connected to any point in A c in the random hypergraph 
T D Del 3 (C). Because of the tile-drawing mechanism, {A va A c } = {Aa«.a c > 0} is also the 
event that there exists a point in £ fl A connected to infinity in T fl Del 3 (£). 

The next Proposition is the key argument why percolation for the random cluster measures 
leads to a break of symmetry in the mark distribution, and as such it is standard and well 
known ( |GH96j . 

Proposition 3.1. For any measurable AcA, 

J {qN A A - N a ) dyv,^ = (q - 1) 


J AA«-a c dCf\ iW . 


Proof. 


This is proved in [Eye 14 


Lemma 2.17] following ideas in G1196 . 


□ 


For purely Delaunay edge potentials an analogous Delaunay random-cluster model exists for 
edge percolation. Then T R 2 is replaced by £*2 = {e = {x, i/}cK 2 :i/|/} and the tile drawing 
mechanism by an edge drawing mechanism (for details see [ BBD04j or [GH96] ). 
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4. Tile-Percolation 


We establish the existence of tile percolation for the Delaunay random-cluster measure 
when z and the parameter /3 are sufficiently large. Note that for any A C l 2 we write 

Na«-oo(C? T) = ff{x G C n A: x belongs to an oo — cluster of (C, T D T^)}. 

The key step in our results is the following percolation result. 


Proposition 4.1. Suppose all the assumptions hold and that z and f3 are sufficiently large. 
Suppose that A is a finite union of cells A k ,i defined in (A.l) in Appendix A. Then there exists 
e > 0 such that 


N, 


A-fAoo 


d C A ,u > £ 


for any cell A = A k>h any finite union A of cells and any admissible pseudo-periodic boundary 
condition u> G 


Proof of Proposition 4.1 We split the proof in several steps and Lemmata below. Our 


strategy to establish percolation in the Delaunay random-cluster model is to compare it to site 
percolation on Z 2 . We adapt here the strategy in [ GH96] to our cases in the following steps. 
First we employ a coarse-graining strategy to relate each site (k, l) G Z 2 to a cell which is 
a union of parallelotopes. In order to establish site percolation we need to define when cells 
are good (open) and when two neighbouring cells are linked once they are open. This link 
establishes then an open connection in our Delaunay hypergraph structures Del m (C),m = 2,3. 

Step 1: Coarse graining. 

Let A = A„ C 


be the parallclotope given as the finite union of cells (A.l) with side length 


f i.e., 


An = (J A kt i and A kji = (J A^, 

(k,l)£{—n,...,n} 2 0 

where A^ are parallelotopes with side length £/9. The central band of two neighbouring cells 
A k: i and A fc+lii is defined by 

ciwu - (U^s* 4 ) u (U^„). 

2—0 2—0 


We now define two properties our Delaunay random-cluster measures need to satisfy in order 
to establish percolation in the Delaunay random-cluster models and therefore non-uniqueness 
of Gibbs measures for certain parameters. Given ( e let H k:k+ i,;(C) be the subset of all 

tiles of Del 3 (C) whose circ nm circles have a nonempty intersection with CB k _k+i,i- Define the 
event F k j that all small cells Aff in A k j contain at least one point of f. i.e., 

Fu =n {#« n Nk > i}. 

i,j =0 

We say the Delaunay random-cluster measure C\ jU} is coarse grain ready (CGR) if there 
are 0 < G < £* < oc such that for all scales £ G (£*,£*) there is a continuous, strictly increasing 
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Figure 1: The central band CBk-.k+i,i is the dark shaded strip (union of 10 smaller cells) 
connecting two neighbouring cells. Some hyperedges are drawn which meet the central band. 


function g of the parameter (3 such that for all neighbouring cells A k,i and A^+i^ and for all 
configurations C £ Fk,i i~l F k+ ij, 

0(a(r)) > g{(3) for all r G H k:k+U (Q. 


To define the second property we need to consider the distribution of the points given by the 
marginal distribution 

,w = x T) 

of the Delaunay random-cluster measure on Note that (3.3) can be written as 

n K(Q,T). /JJA 

C A ,„(dC, dT) = M A ,„(dC)Ai(dT). ,,<’’ A (dT) = 

We define h A to be the Radon-Nikodym density of with respect to P Au , i.e., for ( G Da,u>j 

M0 := Za(C _ 1 / y' (CT W(dr). 


We say the Delaunay random-cluster measure C(y w has bounded Papangelou conditional 
intensity (BPI) if there exists 5 = 5('0,0) > 0 such that for M Aj£J -almost all £ G and 

e A \ C, 


MC U{x 0 }) _ 5 

MO - 9 


(4.1) 


Suppose our Delaunay random-cluster model satisfies both the (CGR)- and the (BPI)- 
property. An important component of our coarse graining method is to estimate the conditional 
probability that at least one point lies inside some A C A. In the following we denote by MyAy 
the conditional distribution (relative to M AtU ) in A given that the configuration in A\ A is equal 
to u / with u' \ A = cu. The details of the construction of the regular conditional probability 
distribution can be found in |Eyel4| . 

Fix e = |(1 -p). site) (Z 2 )). 

Lemma 4.2. Let A as above and let £ > 18Ao- There are constants c = c(5o) and zo = c(6o)q 5 
such that for z > z 0 and for all admissible pseudo-periodic boundary conditions co G 

M w (#{Cnv}>i)>i-|- 
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for all cells V = (k, l ) G {— n ,..., 7i} 2 , i,j = 0,..., 8, and for any admissible boundary 

condition u / G fly with uj' \ A = oj. 


Proof. The statement follows immediately from (BPI) and [E yel4[ IBB D04J using that the 
particle distribution is governed by a hard-core-potential if with range So- Let Vo be the cell 
V reduced by a boundary layer of thickness <5 0 . In particular there is z 0 such that 

,5 


M w (#{(nv} = o)< 


q 


^IVn 


£ 

< 81 


(4.2) 


for all z > Zq. 

£ > 185 0 . 


For the proof the coarse graining length scale £ must be chosen such that 

□ 


Step 2: Random-cluster measure C\ )U . We find a measure C\ )i0 which is stochastically 
smaller than C\ t0J . Then using coarse graining and comparison to site percolation on 1? we 
establish percolation for C\. UJ . Percolation for C\. UJ then implies percolation for the original 
Ca jW . We base the definition of the measure C\ )U on a coarse graining method originally 
introduced in |Hag00| and then used in BBI)() 1 . 

For given ( G let 

Del^C) = {t G Del 3 (C): (j>{a{r)) > g{/3)} 
and let be the distribution of the random set {r G T<y v(t) = 1} with 

l-exp{— g{p)} 


(v(t) = 1 ) = p = 


-1 


Del^(C)( T ) 


(4.3) 


l + (g-l)exp{-^(/3)} 

Note that Jiq depends on A only via the configuration ( G Oa,w Note also that p is increasing 
in /3, although to reduce excessive notation, we don’t explicitly write this. It is easy to show 
that he (see |Eyel4[ Lemma 3.11]). Hence, C\ )U )p C\ )U . As site percolation implies 

hyperedge percolation we consider site percolation given by (4.3), that is, we open vertices in 


Deli(C) independently of each other with probability p. Formally this is defined as follows. We 
let A<^a be the distribution of the random mark vector a G FA where (fT x )x£( are Bernoulli 
random variables satisfying 

P(3y, = 1) =Ffl De i* (c) (a;) 




± !) = 1 -plDelKC)^)* 


where p is given in (4.3) and Del^(£) is the set of points that build the hyperedges of Del^C). 
Then the site percolation process is defined by the measure 

Ci“*>(dC) = M A ,„(dC)A c , A (dC). 

Step 3: Site percolation. 

Given (k,l) G {— n,... , n} 2 we need to estimate the probability of the event that all cells 
A^j, i,j = 0,..., 8, contain at least one point. 


Pick some V = Ay) and choose according to 


Lemma 4.2 z so large that for all admissible configurations u/ G Oy with \ A = a;, we obtain 


M w (#{Cnv} = o)<-. 


Hence, for u>' G we obtain 
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We now establish percolation for the random-cluster measure Ca,w 


Lemma 4.3. There exists /3 q > 0 such that for (3 > /3 0 and, z > zo where zq is given in 
Lemma \f.%\ there is a c> 0 such that 


+■ 


-> A 0 }) > c> 0 


for any A = A k j C A ,(k,l) G {— n, ..., n} 2 , and any pseudo-periodic admissible boundary 
condition cu G 


Proof. Pick ( k , l ) G {— n ,..., n } 2 and consider the event C k) i that each cell A^\ has at least 
one point and all points in A k j fl Del^(C) are carrying mark 1, 

C k j = {C e : C e F kyt and a c (x) = 1 for all x G A k j fl Del^(C}- 


A cell A kt i is declared to be ’good’ or ’open’ if C k j occurs. Each vertex x G Del*(£) is open 
with probability p (see (4.4)). The hard-core background potential if provides an upper bound 
on the number of points in A k ,i for any configuration (. Thus with 

#{A*,i n Del*(C)} < v ~ ' :: u/ =: M 

we obtain 

CZ\C tf )> J M AAh ^,(dQi F J()p M >a-c)p 

Using (CGR) we can now find f3 0 such that for all f3 > f3 0 we have 

p > (1 — e) 


(£ + 25 0 ) 2 
n5l 


-M 


and thus 

Ct“’(Cu) > (1 - e) 2 > 1 - 2 £ = p«“>(Z 2 ). 

The key point is that this property holds for any pair (k, l) independently from any other pair 
of box labels. Standard results for site percolation on Z 2 now ensure that the probability of 
a chain of ’good’ cells A m;r from A k ,i to A c is strictly positive. To see this note that each 
site (k,l) G Z 2 is open once C k ,i happens. Therefore there is strictly positive probability 
connecting some cell to the boundary of A via open sites (’good’ cells). It remains to check 
whether this chain also ensures an open connection (path) from to A c within the Delaunay 
graph Del 3 (C). Our previous def ini tion of the central band is designed in such a way that 
two neighbouring cells which are both open (’good’) guarantee an open path using hyperedges 
which have an intersection with the central band. Suppose £ G C k ,i O Ck+ij, that is, the 
neighbouring cells are ’good’. Since C\ tUJ satisfies (CGR) it follows that </>(a(r)) > g(/3) > 0 for 
all r G H k:k+U ((). Thus H k:k+1> i(() C Delg(C) and henceforth cr c (x) = 1 for all x G H k:k+1> i(C). 
In other words, for all (k,l) G {—n,... ,n} 2 the restriction of Del 3 (C) to H k:k+ i^(C) equals the 
restriction of Del 3 (£) to H k:k+ 1,; (C)■ Hence we can connect to A^ u in the graph Del 3 (0 
inside A k) i U A^+i^ (via central band hyperedges) connecting points of mark 1. Hence, there is 
c > 0 with C\ tU j({Ak y i <— * A c }) > c > 0. □ 


5. Proofs 

This section is devoted to the proof of our results. We first establish the existence of a Gibbs 
measure in our models. 







PHASE TRANSITIONS 


15 


5.1. Existence of Gibbs measures. To show the existence of Gibbs measures (Proposi¬ 


tion 2.2) for all our Delaunay Potts models we will show that they satisfy adapted versions 


of the conditions necessary for the existence outlined in |DDG 12 j . For Delaunay hypergraph 
structures the authors of 1)1)0 12] have given three general conditions hyperedge potentials 
need to satisfy to ensure existence of Gibbs measures. The potentials 0 and 0 of all our models 
depend solely on the Delaunay hyperedges Del m (£), m — 2 or m — 3 , of a marked configu¬ 
ration <0 Every marked hyperedge 77 G Del 2 (C) has the so-called finite horizon B(rj,Q where 
B(r), () is the open ball with 8.8(77, () D ( = 77 that contains no points of (. The finite-horizon 
property of a general hyperedge potential 92: fiy x —» M says that for each pair (77, £) with 
77 G Del m (£) there exists some A d M 2 such that for the pair (77, £) with 77 G Del m (£) we 
have that 92(77, £) = 92(77, £) when ( = ( on A e B( The hyperedge potential 92 in all 
our models is the sum of two potentials, the background potential -0 and the type potential 
0 . As all our potentials depend solely on the single hyperedge they share this finite-horizon 
property. Thus they satisfy the range condition (R) in [ DDG 12 ] . see [DDG 12 , Proposition 4.1 
& 4 . 3 ]. The second requirement is the stability condition. A hyperedge potential 92 is called 
stable if there is a lower bound for the Hamiltonian for any Adi 2 . This stability condition 
(S) is satisfied for all our potentials 0 and 0 as they are all positive. The third condition to be 
checked is a partial complementary upper bound for the Hamiltonian in any A d M 2 . This is 
a bit more involved and we shall first define appropriate configurations, the so-called pseudo- 
periodic marked configurations. We consider the partition of M 2 as given in Appendix |Aj We 
let .8(0, r) be an open ball of radius r < p 0 £ where we choose p 0 > 0 sufficiently small such 
that .8(0, r) C A 0 ,o- Note that 

B r : = {( G Da 0 ,0 : C — {^} for some x G -B(0,r)} 
is a measurable set of Da 0 , 0 \ { 0 }- Then 

T r = {uGfl: 0mz{u A kl ) E B r for all (k,l) G Z 2 } 


is the set of pseudo-periodic configurations (A.2). Note that these configurations are not 


marked. The reason is that when a point is shifted its mark remains unchanged. Thus we 
define the set of pseudo-periodic marked configurations as 


r r = {uj = 






U! 


(i) 


G T r for all i G E}. 


The required control of the Hamiltonian from above will be achieved by the following proper¬ 
ties. As all our hyperedge potentials depend only on the single hyperedge the so-called uniform 
confinement (see }DDGl2j ) is trivially satisfied. In addition we need the uniform summability 
for any hyperedge potential, that is, 

9 9 ( r 7>C) . 

c r ■= sup y — < 00 , 

cer r ^ J #7? 

77eDel m (C): r;nA^0 " 

where rj = {(k,l) G Z 2 : 7 / D A kyi 0 0}. Note that the length l{rj) of any 77 G Del 2 (C) H A 
when Q is any pseudo-periodic configuration satisfies £(1 — 2p 0 ) < £(rj) < 1(1 + 2p 0 ). As all our 
models have the background hard-core potential 0 with range (0 we need to choose the scale l 
for the partition such that £(1 — 2p 0 ) > We compute the constant c r for each model. Our 
chosen partition implies that maximal six Delaunay hyperedges exists. Furthermore, rj = 2 
when 77 G Del 2 (C) and rj = 3 when 77 G Del 3 (C)- 
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Triangle model I: As 0(a(r)) < P for any r G De^((^) we obtain 

c r = 2/3, 

and thus uniform summability holds as long as f3 < oo. The so-called strong rigidity is given 
once z\B(Q,r)\ > e Cr , that is, once 


, 2/3 


Z > 


—: Zq(Sq, /3). 


7Tpl(185 0 ) 2 

Triangle model II: As 0(a(r)) < log(l + /3(tt/3) 3 ) for any r G Del 3 (£) we obtain 

c r = 21og(l + /3 (tt/3) 3 ) 

and thus uniform summability. We obtain strong rigidity once 

z> ^g(iW (l+(/i,r3)/27)2=:z ° (5o ' /?) ' 

Edge model: As <f{£{rj)) < log(l + /3(5 0 /£(l — 2p 0 )) 3 ) we obtain 

P 


c r = 3 log 1 + 


(18(1 - 2p 0 )) 3 

and thus uniform summability. We obtain strong rigidity once 

If. 0 \3 


< OO, 


z > 


1 + 


—: ^o(Aq, /?). 


np 2 0 (185 0 ) 2 V ' (18(1 — 2p 0 )) 3 - 

Using (DDG121 Theorem 3.3] and jDDG121 Corollary 3.4] we obtain all the statements in 
Proposition 222 by noting that for any given z > 0 we can pick po and £ sufficiently large such 
that the above conditions on z are satisfied. □ 


5.2. Percolation of Delaunay random cluster models. We show percolation in all Delau¬ 
nay random-cluster models for sufficiently large parameter z and (3. As outlined in Section [4] 
we need to establish the (CGR)-property and the (BPI)-property. 

Lemma 5.1 (CGR - Triangle Model I). For all £ F k j f~l F k -^-\j and t £ (0, 

(M) e {- ? V • -,ri} 2 , 

0(o(r)) > P, 

where the coarse graining scale £ satisfies £ G (185 0 , m 8 0 ) for some m > 18 and a 0 G (0, ^h_). 


Proof. Suppose that 2<5 0 < £ < m8 0 and a 0 G (0, ^^). Since all of the little cells Ap), i,j = 
0,..., 8 with side length £/9, contain at least one point, we have that any open ball of radius at 
least jh\/7£ and centre xq G CBk-.k+\,i has a non-empty intersection with Therefore, for each 
r G H k: k+i,i(C) the circumscribing circle B(r, () has radius less than j^\fl£. Let r = {a, b, c } be 
such a triangle. Without loss of generality, let a(r) be the angle acb and l be the arc length of 
the arc on dB(r , () between a and b. Let x be the centre of B(t, £). It follows that axb = 2a (r) 
and l = 2 ra(r) where r is the radius of B(r,Q- By the hard-core condition l > \a — b\ > So- 
Combing these facts with r < ~^V7£ gives 


r \ ^ 9h° 9 

a{r) > —^=- n > —f =— > a 0l 


y/7 £ \/7m 

and thus f>(a(r)) > P for all r G H k: k+i,i(()- We obtain (CGR) by choosing g(P) = p. 


□ 
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Lemma 5.2 (CGR - Triangle Model II). For all G F ki i f~l F k - i-i ? z and, t G (0, 

(M) G {-n,... ,n} 2 , 

</>(o(t)) > log(l + ^«o), 

where the coarse graining scale £ satisfies £ G (18<5o,m<5 0 ) f or some m > 18. 


Proof. With at(r) > at 0 from Lemma 5.1 we obtain the statement. Note that the lower 
bound is strictly increasing in the parameter /3. Thus we conclude with g(/3) = log (l+^Og). □ 


Lemma 5.3 (CGR - Edge Model). For all G F k j f~l F k ^-^j aud oil cdcjcs rj G DgI 2 ^^) such 
that r) C t for some r G Hk : k+i,i((), (&, l ) E { —..., n} 2 , 

> log (l + l3/(9V7m) 3 ), 

where the coarse graining scale £ satisfies £ G (185o,m5 0 ) for some m > 18. 


PROOF. From the proof of Lemma |5.1| we know that the radius of the circumcircle of any 
triangle r G Hk : k+ yz(C) I s l ess than j$V7£ and thus any Delaunay edge in Del 2 (C) is no longer 
than \'Jl£. We thus conclude with the lower bound noting that this bound is strictly increasing 
in /3 by choosing g(f$) = log (l + /(9%/Tm) 3 ). □ 


In order to show (BPI) for all our models, we need to investigate the geometry of the 
Delaunay triangulation Del 3 (C), and in particular what happens to it when we augment ( with 
a new point xq f (. Some hyperedges may be destroyed, some are created, and some remain. 
This process is well described in |Lis94j . We give a brief account here for the convenience of 
the reader. We insert the point Xo into one of the triangles r in Del, 3 (<C). We then create three 
new edges that join xq to each of the three vertices of r. This creates three new triangles, 
and destroys one. We now need to verify that the new triangles each satisfy the Delaunay 
condition (2.1), that is, that their circumscribing balls contain no points of (. If this condition 
is satisfied the new triangle remains, if it is not satisfied, then there is a point ij G ( inside 
the circumscribing ball. We remove the edge not connected to x 0 , and replace it by an edge 
connecting xq and X\. This results in the creation of two new triangles. Each of these triangles 
must be checked as above and the process continues. Once all triangles satisfy the Delaunay 
condition, we arrive at the Delaunay triangulation Del 3 (() U {xq}). Let 


= Del 3 (C) n Delate U {^o}), 

T Zi = De| 3 (C u No}) \ Del 3 (C) = Delate U {*„}) \ T}">, (5.1) 

r k< = Dela(e) \ Delate U {*„}) = Delate) \ T<“", 
be the set of exterior, created, and destroyed triangles respectively, see figure [2j 
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Note that any new triangle must contain xo, i.e., 

T io,C = i T e De| 3(C U {x 0 }): r fl x 0 = x 0 }. 

We let ^~ 0 \, l^xoO an d l- l xo l be fl ie tile (triangle) drawing mechanisms on T^, T^, and T^, 
respectively, which are derived from the tile drawing measure /j,^\ in Section [3j We need to 
study the change of KT),T C T ; - C 7 when adding the point x 0 ^ (. Adding a point x 0 to ( 
without considering the change to T will always increase the number of connected components 
by one. On the other hand, the augmentation of a single triangle r to T can result in the 
connection of a maximum of three different connected components, leaving one. Therefore, 

K(C U M,T) - A'(C,T) = 1, 

-2 < /i(C.TUr) — A'(C,T) < 0. 

Lemma 5.4 (BPI - Triangle Model I). For A d M 2 and pseudo-periodic admissible boundary 
conditions u G and M\ jOJ -almost all ( G and xq G ( \ A, 

MC U {x 0 }) 


MO 


> q 


47r 
’“0 , 


PROOF. The tile drawing mechanism /a^.a opens only triangles r with minimal angle a(r) > 
do- We collect all triangles touching A and having minimal angle greater than «o into the set 
T*. Using this fact and Del 3 (<C U {x 0 }) = T^ U T^, we get 

MC U {i„}) St- 9' f(CU( " l ' T Vcu( w) .A( dr ) 


MO 


St* 9 k «' t) hA<IT) 


St*„t>-"> St*„t<« 






W-) Wh <i m ' J (d r 4 )^(dr 3 ) 


^O.c 


^O.c 


Note that the maximal number of triangles in T 2 is Therefore, by (5.2), we conclude with 


K{ C U {x 0 }, T A U T 2 ) - K(C, Ti) > - —, 

do 

AT(C, T 3 U Ti) - K(C, T 3 ) < 0. 


For our triangle model 11 we shall obtain a lower bound for the difference 

iF(CU{x 0 },TiUT)-iF(C,Ti), 

where T\ C and T C T^. Application of (5.2) would result in 

K(( U {x 0 }, T x U T) - K{ C, TO > -2#{T}. 


□ 


(5.3) 

For the triangle model 1 it was easy to bound the number of elements in T via the minimal 
angle condition. In the triangle model II we only obtain a bound on the expectation of that 
number. We obtain the bound on the expectation. 

Lemma 5.5. Let A d M 2 and ( G £Ia,w Then 

[ #{DOdT) < 2x(l + V 2 ) 


for x 0 G A\(. 
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Proof. As any r G touches A we get with the tile drawing probability p(r) in (3.2), 


#{T}<; c (dT)= £ (l_ e -*(«W)). 


t£T 


(+) 


^0,C 


Write Tj ~'. as a union of the disjoint sets 


Hi = R 6 : q(t) > 1}, 

^ = {re ^xoi: jj; < «( r ) < JATtM ^ 2- 


k 


Note that #{-£/*,} < 2nk for all fceN. Thus 


/ #m<; c (dr) = EE f 1 - e ^ l “ ,r)) ) 

OO 7 -j 

S 2tr(l + ^E ( *_ 1) . + 1 ) S ^ + 3^)’ 

fc =2 V ’ 

where we used that 

OO 7 OO ^ o 

y fc < 2 y- = -. 

(k — l) 3 k 2 3 

k =2 v 7 fc=l 

□ 

We then get (BPI) using Jensen’s inequality. 

Lemma 5.6 (BPI - Triangle Model II). For A d M 2 and pseudo-periodic admissible bound¬ 
ary conditions lu G and M A ^-almost all ( G Oa )CJ and xq G C \ A, 

MC U {^0}) > 47r(l+|/37r 2 ) 

MO - 9 


Proof. We proceed as in Lemma 5.4 Jensen’s inequality, (5.3), and Lemma [575] provide the 
lower bound 

/ <(K (Cu(x» ) .T,UT)- J e(C,T 1 ), 1 (+> {(dT) > f iT) 

> ? -2/#{T}^+) c (dT) > g -4^(l+±/^ 2 ). 

□ 

For the edge model we define Delaunay edge sets E 1 ^, and in the same way as done 

for triangles in ( 5. 1| ) (see also [BBD04J). As outlined in Section [3] an edge drawing measure 
is dehned in the same way as the tile drawing measure [ 1 q )A in (3.2), that is, 


Pa(v) = 


(1 - eWM) H D ei 2( c)(h) if V e E R 2 \ A A c, 


llDel 2 (C)(h) 


if v G E A c, 


where £{rj) is the length of the Delaunay edge r/. Then /W ^ denotes the corresponding edge 
drawing measure on E ( +k. We need to bound the expectation of the edges in E^. 
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Lemma 5.7. Let A d M 2 and ( G Oa, w - Then 

f #(£} O dB ) < 4(9 + i/?* 2 ) 

for x 0 G A VC- 

Proof. Due to the hard-core background potential we obtain almost surely with respect to 
Ma,gj the following upper bound on the number of points y in ( with \y — x 0 \ < r, 

\B(x 0 ,r) n Cl < 4 (~T^) 2 > 

and thus for r = 2So, 4J 0 , •.., 

\B(x 0 ,2S 0 ) n Cl < 36 =: &i, 

\B(x 0 ,2nS 0 ) flCI < 4(2n + l) 2 =: b n ,n > 1. 

The number of edges in Ef^^ that have length less than r is bounded by \B{xq, r) fl C |? however, 
as r increases, this bound grows quadratically. On the other hand, according to the edge 
potential </>, edges with large length are less likely to survive the pA-thinning process of E^^ 
than their counterparts with small length. We tradeoff these facts in the following. Let 

OO 

M 2 = A n with A 0 = B(x o, So) and A n = B(x o, 2 nSo) \ A n _i, n > 2. 

n =0 

Now let E be the pA-thinning of Ef^^ and dehne 

E n = {rj = {T 0 ,a;} G E: \x 0 — x\ < 2nS 0 }. 

Clearly, f #{E 1 }is^(dE) < 36. Now for n — 2 at most 36 points are in distance less than 2<5 0 
from Xo, whereas the remaining points lie in the annulus A 2j and due to the fact that Pa is a 
decreasing function of the length (distance), they have at most probability p\(2So) of sharing 
an edge with xo in E. Therefore, 

/ n —1 n —1 

#{£“} <ydB) < 36 + y> +1 - b t )p A (2kS 0 ) < 4(9 + 8 £ kp\(2k8 0 )). 

k=1 k= 1 

Using 

PA (H 0 ) = l^e-«“.) = F T_, 

and the monotone convergence theorem, we finally obtain the bound 

/ #{E} C c (dS) < 4(9 + PY. p) = 49 + -JA' 

J k= 1 

□ 

Having this bound we immediately obtain the (BPI) for the edge model using Jensen’s 
inequality. 

Lemma 5.8 (BPI - Edge Model). For A <<= M 2 and pseudo-periodic admissible boundary 
conditions u G and M^-almost all ( G Oa,u and xo G ( \ A, 

MC U W) ^ —4(9+g/37r 2 ) 

Mo - q 
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We finish the proof of percolation in all Delaunay random-cluster models: 


Triangle Model I: We notice that (BPI) holds for 6 = — (Lemma 5.4h . We now pick 


z 0 = Zo(do,ao) according to Lemma 4.2 and from (CGR) in Lemma 5.1 we obtain /3 0 = /3o(ho) 
such that Lemma |4.3 holds for z > Zo and /3 > /3q. Henceforth, we have that Proposition 4T 
holds and with Section 5.3 we conclude with Theorem 2.6[a). 


Triangle Model II: We see that (BPI) holds for <5 = Att{1 + ^(5tt 2 ) and observe that it depends 
on /3. Thus we first pick /3 0 such that p > (1 — e)^ in the proof of Lemma 


4.3 


We then pick 


zo = zo(6 0 ,j8) according to Lemma 4.2 i.e., (4.2), and from (CGR) in Lemma 5.2 we obtain 


Lemma 4.3 such that Proposition 4.1 is established for all z > zq, and finally Theorem |2.6[b). 


Edge Model: We see that (BPI) holds with 5 = 4(9 + j-/3n 2 ) and observe that it depends 
on /L Thus we proceed as before for the triangle model II. Note that we need an adaptation 


of Proposition 4.1 to the hypergraph structure Del 2 (C)- This follows as outlined in |Eyel4|| | or 
[ BBD04] . We thus have established Theorem 2.4 


5.3. Breaking of the symmetry of the mark distribution. In this section we complete 
the proof of Theorem 2.4 and Theorem 2.6 by analysing the Gibbs distributions 7 a )W hi the limit 
A 'l M 2 . We pick an admissible boundary condition co € Q* Vn and let u> = (u\ A„, 0,..., 0) be 
the admissible monochromatic boundary condition such that l c £ O^. We write for 7 a„,<xj 
for ease of notation and we let P n be the probability measure on fl relative to which the marked 
configurations in distinct parallclotopes A n + (2n + 1 )M(k, /), ( k , l) £ Z 2 , are independent with 
identical distribution 7 n . As we are dealing with a cell structure for the partition of M 2 we 
confine ourself first to lattice shifts when we employ spatial averaging. Thus, 


P n = 


2 n + 1 ^ 

n,...,n } 2 


Pn ° 


By the periodicity of P n the measure P n is Z 2 -shift-invariant. The proof in [ DDG12I Chapter 
5] shows that (P n ) n > i has a subsequence which converges with respect to the topology of 
local convergence to some P £ A4i(f2). As outlined in II )1 )G 12 it is difficult to show that 
P is concentrated on admissible configurations. As P is non-degenerate the proof in |DDG121 
Chapter 5] shows that P = P(-|{0} c ) is a Gibbs measure with P({0}) = 0. In order to obtain 
an M 2 -shift-invariant Gibbs measure one needs to apply another averaging, 


p d) _ 


' Ao,o 


P O 6j x dx. 


Applying Propositions 3.1 and 4.1, we see that for A = A 0 0 , 


{qN A ,i ~ N a ) dP n > 77 —y 

Zn + 1 


E 


N, 


A k i-H-oo cICa^^ 


(fc,0e{— n,...,n} 2 


>(q- l)e- 



N a ) dP (1) > 0, 


Thus 
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and using the symmetry of P (1) we observe the following break of symmetry in the expected 
density of particles of type 1 and of any other different type, i.e., 

Pl (p m ) > P2 (.p m ) = ■ ■ ■ = P ,(p m ), 

where p s (P (1) ) = l/|A|E P (i) [-/Va )S ], s G E. We conclude with our statement as in 1( 11196 by 
showing that the matrix 

CT (0 ))., (6E 

is regular, where P (f) is obtained from P (1) by swapping the role of 1 and t. 

Appendix 

A. Pseudo-periodic configurations 

We define pseudo-periodic configurations as in |DDG12j . We obtain a partition of M 2 which 
is adapted to the Delaunay tessellation. Pick a length scale l > 0 and consider the matrix 



Note that [M/ — £,i — 1,2, and Z(M 1; M 2 ) = n/3. For each ( k,l ) G Z 2 we define the cell 

A M = {Mx G M 2 : X - (k, l) G [-1/2,1/2) 2 }. (A.l) 

These cells together constitute a periodic partition of M 2 into parallelotopes. Let B be a 
measurable set of Da 0 , 0 \ I 0 } and 

r = {uefl: e Mz (u Aktl ) G B for all (k, l ) G Z 2 } (A.2) 

the set of all configurations whose restriction to an arbitrary cell, when shifted back to Ao,o, 
belongs to B. Elements of T are called pseudo-periodic configurations. We define marked 
pseudo-periodic configurations in an analogous way. 

B. Topology of local convergence 

We write Adf(O) (resp. A / lf(D)) for the set of all shift-invariant probability measures on 
(Q, J-) (resp. (0,P)). A measurable function /: Q —> M is called local and tame if 

/M = f(w A ) and |/M| < aN A (u) + b 

for all uj G fl and some Ail 2 and suitable constants a, b > 0. Let L be the set of all local and 
tame functions. The topology of local convergence, or Z-topology, on A4f (fi) is then defined 
as the weak* topology induced by A, i.e., as the smallest topology for which the mappings 
Pig / fdP with / G A are continuous. 
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